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We obtained steady solutions to the two-dimensional Boussinesq approximation equations without 
mean temperature gradient. This system is referred to as free convection in this paper. Under an 
external flow described by the stream function \P = —Ayf(x), two types of steady solutions are 
found depending on the boundary conditions. One is kept steady by the balance between the strain 
of \P and the diffusion. The solution is similar to the Burgers vortex layer solution. The other is done 
by the balance between vorticity induced by the buoyancy and vorticity flux caused by the external 
flow. Detailed argument on these two balances is presented for f(x) = x. Then two examples other 
than f(x) — x are shown to have either of the two balancing mechanism. We discuss the relation 
between these solutions and long-lived fine scale coherent structures observed in direct numerical 
simulations of two-dimensional free convection turbulence. 

PACS numbers: 47.27.-i, 47.27.Te, 47.32. Cc 



I. INTRODUCTION 

Rayleigh-Benard convection is very rich in nature. It 
has been a source of inspiration in studies of universality 
in the route to turbulence B. At the same time, there 
have been extensive researches both experimental and 
theoretical on turbulent states of convection [^J. 

Recently, convection turbulence has attracted atten- 
tion as a good example of active scalar systems. One 
interesting observation is scaling behavior of structure 
functions of temperature, active scalar. In this con- 
text, a two-dimensional Boussinesq convection model has 
been studied by direct numerical simulations (DNS) in 
a doubly periodic square domain. The results indicate 
that the scaling exponents of the temperature structure 
functions saturate as the order goes to infinity & Q. 
It has been pointed out that coherent structures like 
plumes or sharp interfaces (shocks) between hot and 
cold regions can play an important role in this behav- 
ior. Moreover it is shown phenomenologically |^| and 
numerically |6) that the temperature variance (entropy 
||) cascades from larger scales to smaller scales in the 
turbulent two-dimensional Boussinesq system (referred 
to as two-dimensional free convection, 2DFC) when the 
turbulent state is maintained by a large scale tempera- 
ture forcing. However, the relation between the entropy 
cascade and the dynamics of coherent structures, if any 
at all, is still far from understood. 

In the three-dimensional Navier-Stokes (3DNS) turbu- 
lence, vortical fine coherent structures, so-called worms, 
are observed. These worms are well known to be approx- 
imated locally by the Burgers vortex tube. In this sense, 
2DFC turbulence is quite similar to 3DNS turbulence, 
although the former is more feasible than the latter. We 
expect that the understanding of the former may help 



us understand the latter and also elucidate the universal 
characteristics of turbulence. 

The motivation of this paper is to explain some proper- 
ties of the long-lived coherent structures observed in DNS 
in terms of a class of steady solutions of 2DFC. For this 
we search for steady solutions to the incompressible two- 
dimensional Boussinesq approximation equations with- 
out mean temperature gradient (2DFC): 

d t T+{u-V)T = kAT, (1) 
Vn 

d t u + (u ■ V) u = + agTe + vAu, (2) 

Po 

V-it = 0. (3) 

Here k, po, a, g, and v are the molecular diffusivity, the 
mean density of the fluid (we take po to be unity for 
simplicity), the thermal expansion coefficient, the gravi- 
tational acceleration and the kinematic viscosity, respec- 
tively. The vector e is the unit vector in the direction 
opposite to the gravity. In DNS of free convection turbu- 
lence, to keep the system statistically stationary, a large- 
scale forcing term is added to Eq. (|l|) , but the details are 
not mentioned here. 

In the turbulent temperature field, long-lived sharp in- 
terfaces between hot and cold regions, i.e. shock fronts 
are often formed and regarded as coherent structures. 
Thus coherent structures are highly elongated; The typi- 
cal length is the order of the integral scale and the typical 
width the order of ten times of Kolmogorov dissipation 
length scale. We expect that steady solutions to Eqs. (|I|)- 
(^) have common properties with the coherent structures 
formed spontaneously in DNS. 

To deal with temperature shocks more directly, we use 
the following vector quantity [Q, || : 

X = (d y T, -d x T) . (4) 

We call x T-vorticity. T-vorticity obeys the following 
evolution equation similar to the three-dimensional vor- 
ticity equation: 
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dtx + (« ' V) x = (X ■ V) u + kA X - (5) 
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Thus T-vorticity is expected to play an essential role 
in turbulence like vorticity in 3DNS. In contrast with 
T-vorticity, the governing equation for the vorticity 
u)(x, y,t) in 2DFC can be written 

d t Lu + (u-V)lo = ag{V x (Te)} z + i/Aw. (6) 

Here {-} z denotes taking the z-component. Note that 
because of the buoyancy term, i.e. the first term in the 
r.h.s. of Eq.(^), vorticity is not a conservative quantity 
even for v = 0. For simplicity, we assume v = k , i.e. 
Prandtl number Pr = v j k = 1 here. 

In the following sections, we look for steady solutions 
to Eqs. (||) and (^) , where the flow is decomposed into the 
external part given by a stream function ^ = —Ayf(x) 
and a response to it. In section II, we deal with the 
case f(x) — x, i.e. a stagnation flow. In Section III, ex- 
tended stagnation flows, f(x) ^ x, are considered, where 
we show two examples of f(x). The realizability of the 
extended external flow is discussed. Concluding remarks 
are made in Section [V. 



II. 



BURGERS VORTEX LAYER TYPE 
SOLUTIONS 



We assume that the system @ and @ is exposed to 
a stagnation flow 



^(x, y) = -Axy, 



(7) 



where A is a positive constant so that we take the direc- 
tions of contraction and expansion as s-axis and y-axis, 
respectively. The angle between ir-axis and e is also a 
parameter and denoted by ip as shown in Fig. [l]. Then 
the unit vector e in Eq. ([j]) is given in this coordinate 
as e = (cos ip, simp). We use this coordinate throughout 
this paper. We further assume the following forms of the 
temperature and velocity fields: 

T(x,y,t) = 6(x), (8) 
u{x,y,t) = (d y W(x,y),-d x W(x,y)) + (o,v(x)) 

= (-Ax, Ay + v(x)). (9) 

In terms of 9 and v, T-vorticity and vorticity are ex- 
pressed as 



u(x,y,t) = oj(x) 



(0,-d x 9(x)) = (0, Xy (x)), (10) 

Mx)__ (ii) 

ax 



Consequently, Eqs. (|^) and (||) are reduced to the fol- 
lowing ordinary differential equations: 



_d_ ( K_dXy 
dx \A dx +XXv 



v d 2 LO 

A dx 2 



dx 



0, 



ag sin ip 
~A ' 



(12) 
(13) 




FIG. 1: The coordinate system fixed on the external stagna- 
tion flow. The direction of the gravity vector g is expressed 
via the angle p. 



After integrating once, Eq.(|l2|) can be written 

XXy = C. 



K_dXy 

A dx 



(14) 



Here C is an integration constant. 

Let us consider the general boundary conditions at the 
origin: 

(15) 
(16) 
(17) 
(18) 

where x' y ( x ) an d u)'(x) denote the derivatives of T- 
vorticity and vorticity respectively; xo^xb^o an d u' 
are constants. These conditions yield the relation C = 
Kx' /A. With n — v the solutions to Eqs.([l2"|) and (|l^ ) 
are obtained 



Xy(x 


= 0) 


= Xo, 


Xy{% 


= 0) 


= Xo, 


uj{x 


= 0) 


= uo, 


Uj'(x 


= 0) 


= ^ 



Xy(x) 
w(x) 



Xo e 



w 



Xo e 



, ag sin p 



* dt, 
d(e~£< 



e^dt 



agsmip 
Xo - A (1 - e 



A 

C rfJ- 

d/j d£ e 
li Jo 



(19) 



Af 2 

2n ? 



(20) 



Here only the last two terms of Eq. (g0|) are the buoyancy 
generated vorticity. The first two are homogeneous solu- 
tions, i.e., the solutions to the vorticity equation ( |l3| ) set- 
ting the r.h.s. (the buoyancy force) zero. When x'o = 0, 
the solution ( jl9| ) which is the same as the Burgers vortex 
layer, takes finite value of temperature at x — ±oo. We 
call this solution finite temperature (FT) type solution. 
In contrast, when x'o 0, the temperature correspond- 
ing to the solution (Q9|) diverges as x — > ±oo. We then 
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FIG. 2: Finite temperature (FT) solutions: the T-vorticity 
Eq.(E3}, the corresponding temperature 8(x) and the vorticity 
Eq. (p4). Here we set the parameters as 0o = 1, A = 1, k = 
1, ag = 1 and 95 = tt/4. 



call this solution infinite temperature (IT) type solution. 
We treat the two cases separately in the following sub- 
sections. 



A. Finite temperature type solutions 



Let us consider the following boundary conditions: 



Xy( x ) 

L0(x) 



(x 
(x 



±00) 
±00) 



(21) 
(22) 



These conditions (|l]) and (22) are equivalent to x'o — 



and lo'q = 0. The general solutions (|19|) and (|20|) are then 
reduced to Burgers vortex layer type solution: 



X y (x) = 9 \l-^-exp (-— x' 



uj(x) 



M exp (_i^ 

A V 27TK V 2K 



ag sm 
A 



Xy { x )> 



(23) 



(24) 



where <9 Q = J^XyiZit = o) d£ = 0(-oo,o) - 0(oo,o). 
In the Burgers vortex layer, this temperature difference 
corresponds to the circulation of the vortex layer. In 
Fig. I we plot the solutions; T-vorticity (p3]), and the 
corresponding temperature 0(x) and vorticity (|24|). As 
mentioned before, the temperature takes finite value at 
x = ±00. Since the temperature 9{x) is a monotonic 
function, Xy( x ) contains a single bump. We call this 
solution Burgers T-vortex layer solution. 

Now we make several remarks on the FT solutions fl23| ) 
and (pi|). First, we note that u)(x) is proportional to 
Xy{ x ) an d sin tp. That is, Xy{ x ) 1S independent of the 
direction of the gravity, but oj(x) is not. This is be- 
cause the former is maintained by the balance between 



the imposed strain and the diffusion of T-vorticity like 
the Burgers vortex layer solution to the 3DNS equations 
and the latter does by the balance between the diffusion 
of vorticity and the buoyancy induced by the temperature 
difference across the Burgers T-vortex layer. Therefore 
vorticity is maximized when the direction of the T-vortex 
layer coincides with that of the gravity. 

Secondly, it is possible to estimate a characteristic time 
of the relaxation to the steady state (|2^) after the tem- 
perature difference <9 is fixed. Under the stagnation flow 
(H), the time-dependent solution to Eq. (||) is obtained 



Xy(x,t) = 4ttk 



1 - e 



-2 At 



2A 



-1/2 



x / Xj/(C,* = 0)exp 



(x- e~ At 2 
An 



M / A A 

&o\ exp x' 

27TK V 2K 



2A 

as t — > 00. (25) 



The characteristic time for the relaxation is thus esti- 
mated as 1/A. 

Thirdly, the characteristic width of the bump of T- 
vorticity and vorticity is A. This agrees with an 

observation about structures of a passive scalar advected 
by a two-dimensional synthetic random velocity field with 
a finite correlation time in Ref. (Of course, vortic- 

ity is not considered in that case.) Although a passive 
scalar with non-zero mean gradient was considered, the 
solution (^3|) may describe locally coherent structures of 
the passive scalar because of the following fe ature s: The 
width of those was found to scale well with \J k/s, where 
s was a root mean square of the rate of strain calculated 
from the low-path filtered velocity. They also showed a 
stagnation flow accompanied by a coherent structure in 
Fig. 7 of Ref. @. 

Finally, concerning the relevance of the FT solutions 
( p3| ) and (^4|) to the coherent structures observed in 
DNS of free convection turbulence, we refer the reader 
to Ref. [fi"l| . Here suffice it to say that the solutions (|23|) 
and ( p4[ ) can give a good local approximation of the co- 
herent structures in DNS, where we compared A with 
the (positive) eigen-value of the velocity-gradient tensor 
although the latter is defined only locally. 

However, it is well-known that, in the limit of k — > 0, 
the temperature dissipation rate per unit length along y- 
axis, ee, of the Burgers (FT type) solution ( p3| ) vanishes. 

That is, ee = J^L, K Xy( x ) 2 dx cx V kA — > as K — * 0. 
It suggests that this Burgers layer type solution is irrele- 
vant in this limit. In DNS, however, the situation can be 
different because we have to consider feedback to A (the 
eigen-value of the velocity-gradient tensor) from other in- 
gredients of the flow. Indeed, the probability distribution 
functions of A(x,y) obtained in DNS of various k's are 
observed to fall into a unique curve under a suitable nor- 
malization. In other words, A(x,y) field may have some 
nontrivial dependency on k. Thus it is conjectured that 
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the dissipation rate of the FT type solution (|23|), V kA, 
would take non-zero value in the limit of vanishing diffu- 
sivity in real flows. Therefore the relevance of the Burgers 
vortex layer type solutions might not be ruled out. The 
detail of the observation on A(x, y) field in DNS will be 
reported elsewhere. 



B. Infinite temperature type solutions 

In this subsection, we consider the case where the tem- 
perature slowly diverges as x — > ±00, which we call infi- 
nite temperature type solutions. In terms of the bound- 
ary conditions at the origin, it corresponds to 




Xy{x = 0) 


= 0, 


(26) 


X U (X = 0) 


= Xo, 


(27) 


oj(x = 0) 


= 0, 


(28) 


J{x = 0) 


= 0, 


(29) 



where xd 7^ 0- We set xo = to examine an IT type 
solution solely. The general solutions ( JlSj ) and (^0|) then 
take the form: 



Xy(x) 

Lo(x) 



1 ^i_- r ^ 

X e 2 * x 



(30) 







ag sin ip 
~A 



Jo Jo Jo 



The boundary conditions for the vorticity (28) and 
are chosen to involve only the buoyancy effect on vor- 
ticity since we focus on the vortical structure induced 
by the buoyancy. It should be noted that the homoge- 
neous solutions of Eq. (|13|) a re neglected. We plot the IT 
solutions :T- vorticity Eq.(|30|) ,the corresponding temper- 
ature 9(x) and vorticity Eq.(|3l|) in Fig. [|. 

It should be noted that unlike the Burgers type solu- 
tions T-vorticity and vorticity of IT type solutions take 
the same sign around x _ bumps. This behavior seems to 
contradict physical intuition that hot fluid rises and then 
vorticity with the sign opposite to that of T-vorticity is 
induced. However, vorticity is maintained by the bal- 
ance between the buoyancy and not only the diffusion of 
vorticity but the vorticity transfer by a stagnation flow. 
In fact, if vorticity takes non-zero value with the same 
sign of T-vorticity, say positive, at infinity, vorticity flux 
Axlo through a region with the width S centered at some 
point ,e.g. xq, always increases vorticity in the region. 
The substantial increase in vorticity Q u is given by 



Qui — F u Fd 

a( 5 \ ( 6 

= Alx + -juj{xo + - 

= ASlu(xq) + A5xquj'(xo), 



FIG. 3: Infinite temperature (IT) solutions: the T-vorticity 
Eq.(3C),the corresponding temperature 8(x) and the vorticity 
Eq.(31). Here we set the parameters as A = 1,k = l,ag = 
l,<p = 71-/ 4 . Xo = 1,0(3; = 0) = 0,u' o = and k = 1. The 
dash-dotted line in the vorticity figure is a localized vorticity 
solution Eq.(]33l) obtained by setting uj' = —0.23. 




FIG. 4: Schematic view of vorticity balance in a region with 
the width S centered at xo- F u = A(xq + 8/2)uj(xo + 5/2) and 
Fa = — A(xo — 6/2)u>(xo — 5/2) are vorticity fluxes in x by a 
stagnation flow. F x = 5x y {xo) 019 B ™ v is vorticity induced by 
the buoyancy. 



where the first term in the r.h.s balances with the vor- 
ticity flux in y and the second term does with F Y = 



shown in Fig. |j. That is, vorticity induced by the buoy- 
ancy cancels out the excess of the vorticity flux. On the 
other hand, if the sign of vorticity around infinity is op- 
pose to that of T-vorticity, the excess of the vorticity flux 
cannot be canceled out by the vorticity induced by the 
buoyancy so that the vorticity is not kept steady. 

When vorticity is localized around the origin, i.e, van- 
ishes rapidly as x goes to infinity, the vorticity transfer 
is turned off and the vorticity induced by the buoyancy 
works as the vorticity input which balances with the dif- 
fusion of vorticity. This means that to keep the vorticity 
(32)steady, cj(oo) should be null or take the same sign as 
that of the T-vorticity. Even in the IT case, using the 



homogeneous solution to Eq.(h3 
ized around the origin. In Fig 



vorticity can be local- 
3| with dash-dotted line, 



5x y (xo)ag simp/ A if u>"(xo) is negligible in Eq.(|l3|) as we show this localized vorticity, which consists of the 
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buoyancy-induced vorticity ( [Ii"| ) and the homogeneous 
solution, 



u>(x) = ui'q e «^ <i£ 



Here the value for the parameter ljq is chosen to can- 
cel out the uniform vorticity of the solution (|3l"|). The 
form of this localized response vorticity matches the in- 
tuition that hot fluid goes up and cold fluid down. Such 
a localized vorticity structure associated with T-vorticity 
structure is one of the striking features of active scalar 
turbulence. 

For any non- vanishing value of x'oi Xy( x ) is an °dd 
function having a double bump and consequently the 
temperature 8(x) = 6{x = 0) — f Q Xy{£,)d£, forms a ridge. 
Such temperature ridges are often observed in DNS. How- 
ever, far away from the origin, temperature diverges log- 
arithmically. |Q. We think that the divergence at in- 
finite distances is not necessarily pathological, since we 
are interested only in the local behavior of the structures 
(around the origin). Furthermore, this divergence pre- 
vents the strain of the external flow from squashing the 
temperature ridge. Like the FT type, the typical width 
of the vorticity structure is estimated as y n/A. It is 
noted that this IT solution ( |30| ) is also the Burgers vortex 
layer solution to the three-dimensional barotoropic vor- 
ticity equation. (Recall that the governing equation (|^) 
for x is the same as that of three-dimensional barotropic 
vorticity.) 



III. EXTENSION 



In this section, we consider the following external flow: 



V(x,y) = -Ayf(x), 



(34) 



where f(x) is a function satisfying a condition mentioned 
later. When f(x) = x, this flow becomes a stagnation 
flow. A major difference from the f(x) = x cases is that 
the external flow (S3) has in general non-zero vorticity, 



Q(x,y) = -A¥ = Ayf"(x) 



(35) 



Thus f2 should satisfy the vorticity equation. We will 
discuss the realizability of (HJ) in detail later in this sec- 
tion. 

Under the extended stagnation flow (|34|), we further 
assume functional forms of temperature and velocity as 
follows: 

T{x,y,t) = 6(x), (36) 
u{x,y,t) = {dy&,-d x &) + (o,v{x)) 

= (-Af(x),Ayf'(x)+v(x)). (37) 



Then x an d the total vorticity u> are given by 

X (x) = (0,-e'(x)) = (0, Xy (x)), (38) 
u(x,y) = -A<P + v'(x) = Q{x,y) +cj(x). (39) 

We call uj response vorticity. The velocity field Eq. ( |37| ) 
is uniquely decomposed into the y-dependent and the y- 
independent parts. 

Imposing the general boundary conditions 



X y (x = 0) = xo, 

X' y (x = 0) = x'o, 

uj(x — 0) = c^o, 

uj'(x — 0) = lu'q, 



(40) 
(41) 
(42) 
(43) 



steady solutions to Eqs.(|J) and (^) are easily obtained 
for any f(x) as follows: 



Xy(x) 



Xoe 



+Xo e " J o / e" Ja y ' dt; 



(44) 



Cj(x) = UJq + UJ' Q I I 

lo 



.a r< 



So fWdX 



+ / d(e »Jo J( ' I Xy(V>)e" Jo dfi 



Jo 
3A f e -4/ C /W«« 



lu + uj q I e 
lo 



d( 



+Xo — j — / Ce " J dC 



,agsm(p 
+Xo 



A 



d(e 



dJ d ^s:m^ dil 



(45) 



The last two terms of Eq. (|45|) represent the vorticity in- 
duced by the buoyancy. 

Now we discuss the realizability of the external flow 
([m|) . If we claim that Eq. ( [jo]) is a steady solution of the 
barotropic two-dimensional vorticity equation 



(U ■ V)J? = vAf2, 



(46) 



where U = (d y <P, —dx^). Then the equation for f(x) is 



-f^(x) + f(x)f^(x) - f'(x)f"(x) = 0, 



A 



(47) 



where f^'(x) denotes the j-th order derivative. Equa- 
tion (E7p has been studied by many authors as a model 
of a flow near a rigid wall In Ref. jlB) the au- 

thors, who were interested in the finite time blow-up 
of the solution, dealt with the time-dependent stream 
function — —Ayf(x,t) in a bounded domain and dis- 
cussed asymptotic behavior (t — > oo) of the solution to 
the equation for f(x, t), the Proudman-Johnson equation 
[Q. They showed that every solution decays to zero as 
t — > oo with homogeneous boundary conditions. 
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Here instead of solving Eq. ( f47|), we examine the re- 
adability of the external field fl34|) . We do not think 
that the flow (|34T) s hould be an exact solution of the vor- 
ticity equation (}46D. Rather, we assume that the flow 
( |34| ) is formed by a larger scale flow in the turbulence 
field. That is, from a practical point of view, a steady 
solution subjected to an external flow is no more than a 
local model of coherent structures in the turbulent field, 
so that the external flow need not satisfy Eq.([l6|). Thus 
it may be possible to add an forcing term describing an 
effect of the larger scale motions to Eq.(|46|). Denoting 
the forcing term by F(x, y), the equation for fl (or f(x)) 
can be rewritten 

(U = vAQ + F. (48) 

Then the steady-state equation for the total vorticity j39| ) 
can be written 

(w- V)(f2 + c3) = Xy{x)ag sinip + ^A(Q + u)) + F, (49) 

which yields the equation for the response vorticity Q 



v d 2 Cj 
A dx 2 



dx 



Xy{x) 



ag sin (p 



(50) 



from which the solution ( |45|) is obtained. 

By substituting f(x) — x, it is easily checked that the 
solutions ( [44| ) and ( |45| ) contain Burgers vortex layer type 
solutions. 

In the following, we show two examples of f(x). Each 
has both FT and IT type solutions. As the first example, 
we deal with an external flow 



f(x) = Ax 3 - 2x, 



(51) 



which itself has vorticity £l(x,y) — 24Axy. The stream 
lines of the external flow is shown in Fig. |^. Here we 
assume an additional forcing F{x,y) that keeps the ex- 
ternal field ( |5l| ) steady. ( Equation (^) prescribes the 
forcing as F(x, y) = 192A 2 x 3 y) 

For FT type solutions, we impose the boundary con- 
ditions 



Xy( x 


= 0) 


= Xo 


¥= 0, 


(52) 




= 0) 


= Xo 


= 0, 


(53) 


£j(x 


= 0) 


= ^0 


- 0, 


(54) 


uj(x - 


= 0)' 


= ^0 


= 0. 


(55) 



The conditions for the vorticity Eqs.(|54|) and (|5^) are 
chosen again to reflect only the effect of the buoyancy. 
The solutions (ffih and mn) are then reduced to 



Xy(x) 
uj(x) 



Xoe " (x 
ag sin ip 



Xo- 



A 



(56) 
(57) 



^ 




-1 -lW2 



lW2 1 



FIG. 5: Stream lines of an external flow ^(x) = — y(^x 3 — 
2x). The arrows denote the directions of the flow. Stagnation 
points are {-l/y/2, 0), (0, 0) and (l/v^O). 




FIG. 6: FT type solutions for W(x) = -Ay(4x 3 - 2x): the 
T- vorticity (|56|), the corresponding temperature 8(x) and the 
response vorticity (^). Here we set the parameters as A = 
1, k = 1, ag = 1, (p = 7r/4, xo = 1 and 8(x = 0) = 0. 



For IT type solutions, we use the boundary conditions 

(58) 
(59) 



X y (x = 0) = xo = 0, 

x' y (x = 0)=x'o ¥= 0, 

uj(x = 0) = ujq = 0, 

u>(x = 0) = Q' = 0. 



These solutions are plotted in Fig. 



Hence the solutions ( |44| ) and ((45|) is rewritten 

Jo 

,-4(c 4 -c 2 ) 



(60) 
(61) 



(62) 



Q(x) 



Xo- 



ag sin ip 
~A 



'0 

d( e 



0.25 



© 



-0.25 




FIG. 7: IT type solutions for $>{x) 



-Ay(/[x 3 — 2x): the 
T-vorticity Eg. (|62|) , the corresponding temperature 9(x) and 
the response vorticity Eq.(^3|). Here we set the parameters 
as A = 1, k = 1, ag = 1, <p = 7r/4, x'o — — 0) = 0, and 

uj'o — 0. The temperature changes in a non-monotonic way. 
The dash-dotted line in the Q figure is a localized response 
vorticity obtained by setting u>' — —0.172. 



C 



d£ / d/j, e 



Km 4 



(63) 



We plot the T-vorticity fl62|), the corresponding temper- 
ature 8(x) and the response vorticity ( |63| ) in Fig. |?]. T- 
vorticity (|6l) has a double-bump for any non- vanishing 



x'o- For this shape of the response vorticity (|63|), the 
vorticity-flux balance argument given in the previous sec- 
tion can be also applied. In addition, we can obtain a 
localized u> solution by adding the homogeneous solution 
of Eq.(|50|) in a similar way to the Burgers layer type so- 
lutions. 

As an another example, let us consider the periodic 
flow: 



f(x) = sinx. 



(64) 



This flow has vorticity fl = —Ay sinx. The correspond- 
ing forcing F(x,y) is then v Ay sinx. The streamlines 
of this external flow Eq.(|64|) are shown in Fig. ||, where 
stagnation points are located periodically. 

For FT t ype solutions, we impose the boundary condi- 
tions (|52|)-(|55|). Then the general solutions ( fS| ) and j45| ) 
is expressed as 



Xy(x) 
uj(x) 



ag sin (p 



i x) 



Xo- 



A 



-4(l-cos{) 



(65) 
(66) 



They are plotted in Fig. Multiple bumps in Xy( x ) 
are seen. From the monotonic behavior of the temper- 
ature which is similar to the Burgers T-vortex layer so- 
lution, these solutions ( |65|) and (j6(]) are classified to FT 
solutions. The remarkable feature of this solution is the 



FIG. 8: Stream lines of an external flow ^(x) = —ysinx. 
The arrows denote the directions of the flow. Stagnation 
points are (nn, 0), n = 0, ±1, ±2, •• •. 
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FIG. 9: FT type solution for <P(x) = -Ay sinx: the T- 
vorticity Eq.(p5[), the corresponding temperature 8(x) and 
the response vorticity Eq.(p6|). Here we set the parameters 
as A = 1, k = 1, ag = 1, (p = vr/4, xo = 1, 0(x = 0) = 



staircase-like behavior of 9(x). This kind of staircase- 
like change of the temperature is sometimes observed in 
DNS. Hence the solution ( |65| ) and ( |6^ ) can be useful to 
describe the small wavy behavior embedded on a tem- 
perature front. 

To obtain IT type solutions, we use the boundary con- 
ditions (H|)- 10]). The solutions and are then 
expressed as 



Cj{x) 



Xo e 

ag sin ip 



/ e^ 1 " 008 ^, (67) 



.4 



d(e- 



t(l-cosC) 



x j dill d^ 1 -™^. 
Jo Jo 



(68) 



These IT type solutions are shown in Fig. [H]. Its temper- 
ature behavior is similar to the Burgers IT solution. The 
vorticity-flux balance argument can explain the shape of 
u> in this case as well. However the response vorticity 
( |68| ) does not seem to converge as an ±oo, it would be 
no use considering the localized response vorticity in this 
case. 
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FIG. 10: IT type solutions for W(x) = -Ay sin x; the T- 
vorticity Eg. (pi]) , the corresponding temperature 8(x) and the 
vorticity Eq.fl68|). Here we set the parameters as A — 1, « = 
l,ag = l,ip = tt/4, Xo = 0.1 and 0(a: = 0) = 0. 



IV. CONCLUDING REMARKS 

We have examined a class of steady solutions to the 
two-dimensional Boussinesq approximation equations. 
We believe that these solutions well describe coherent 
structures observed in DNS of 2DFC. Under an extended 
stagnation flow = —Ayf(x), we obtained two types of 
exact solutions; FT and IT types. In both cases, T-vortex 
layers are maintained steady by the balance between the 
strain of the external stagnation flow and the diffusion 
of T-vorticity. Furthermore, in the FT case a pair of T- 
vortex layers with opposite signs are annihilated by the 
compression of the extended stagnation flow without the 
continuous supply of T-vorticity which is maintained by 
the infinity of temperature at x = ±oo. 

For the former (FT) type including the Burgers vortex 
layer type solution when f(x) = x, temperature remains 
finite as x goes to ±oo and both T-vorticity and vorticity 
are localized around the origin. In contrast, for the lat- 
ter (IT) type, the absolute value of temperature should 
diverge as x goes to ±oo as mentioned above; T-vorticity 



tends to zero quite slowly and vorticity is saturated to a 
finite value. If f(x) = x, T-vorticity takes at most two 
bumps. On the other hand, if f(x) ^ x then it may take 
as many bumps as possible, although we only dealt with 
f(x) = Ax 3 — 2x and f(x) = simp where two and infi- 
nite number of bumps are observed, respectively. These 
bumps seem to correspond to fine T-vortex structures 
observed in DNS of 2DFC turbulence. 

So far we have taken for granted the independence of an 
external flow and a response. This view may be reason- 
able if we focus on local behavior of the coherent struc- 
tures. Indeed the local shape of the structures are quite 
similar to the steady solutions. Thus it is suggested that 
the characteristic time scale of the external strain field 
are well separated from that of structure's motions. The 
role of the coherent structures on statistical properties 
of turbulence expected to be universal, is still an open 
question. Since the coherent structures have quite long 
correlation length comparable to the scale of the energy 
(temperature variance, i.e. entropy) containing range, 
their existence may break the locality hypothesis signifi- 
cantly. However, the authors have examined the relative 
diffusion in 2DFC turbulence recently and found that 
the coherent structures rather play an essential role in 
keeping the locality. These contradicting characteristics 
of the coherent structures make research on turbulence 
challenging. 

In this sense, free convection turbulence is quite similar 
to 3DNS turbulence. We believe that coherent struc- 
tures are universal and essential ingredients of turbu- 
lence. Therefore extensive researches on coherent struc- 
tures are required. 
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